Revisiting interferences for measuring and optimizing optical nonlinearities 
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A method based on optical interferences for measuring optical nonlinearities is presented. In a 
proof-of-principle experiment, the technique is applied to the experimental determination of the 
intensity dependence of the photoionization process. It is shown that it can also be used to control 
and optimize the nonlinear process itself at constant input energy. The presented strategy leads to 
enhancements that can reach several orders of magnitude for highly nonlinear processes. 

PACS numbers: 42.65.-k,42.25.Hz,52.50.-b 



I. INTRODUCTION 



Since the discovery of the laser, technological pro- 
gresses have allowed to produce optical sources of con- 
stantly growing power. Laser intensities can now be so 
high that atoms and molecules, when exposed to such 
electric fields, can absorb simultaneously a large number 
of photons, leading to highly nonlinear dynamics and to 
the observation of phenomena such as strong field ion- 
ization or high-order harmonic [H and attosecond pulse 
Q generation. More particularly, at moderate intensity, 
the efficiency of these processes scales as a fixed power of 
the intensity. For instance, the rate of multiphoton ion- 
ization is well approximated by where / is the laser 
intensity and is the minimal number of photons that 
the electronic wavepacket needs to absorb for escaping 
from the ionic attractive potential. As the intensity rises 
(typically higher than a few 10 TW/cm^), the electronic 
structure of the system is strongly distorted allowing the 
occurrence of processes such as tunnel ionization, above- 
threshold ionization, and ionization channel closure. As 
a consequence, the very simple relation linking laser in- 
tensity and ionization does not hold anymore. Instead, 
only an effective nonlinearity, which can also depend on 
the intensity, can be extracted from experiments. For in- 
stance, an effective nonlinearity A^eff —7.5 d, Q has been 
determined for the ionization rate of argon at 800 nm for 
/ :^ 50 TW/cm^, while a nonlinearity of 11 is expected in 
the multiphoton regime. The nonlinearity of a process V 
is generally determined by fitting the measured process 
yield by a function proportional to /^^ff. Unfortunately, 
due to the highly nonlinear nature of the process, the dy- 
namic range needed for an accurate measurement of such 
an effective nonlinearity is hardly achievable experimen- 
tally and leads to large uncertainties in the determina- 
tion of A^eff- In this paper, a simple method to determine 
the nonlinearity of any optical processes is implemented. 
This method is based on optical interferences induced by 
two crossed laser beams and relies only on two measure- 
ments in contrast with standard methods which need to 
reconstruct the dependence of the process as a function of 
the intensity on several decades. The presented method, 
which can be applied to any nonlinear process, is demon- 
strated analytically, numerically, and experimentally in 



the case of ionization. The proposed strategy also pro- 
vides a control and enhance the production efficiency of 
the nonlinear process. In particular, we show that the 
yield can be increased by several order of magnitude and 
that the overall energy needed for the production of any 
nonlinear process in the two-beam geometry is reduced 
up to a factor 2 as compared to the single beam geom- 
etry. In the first section, the principle of the method is 
presented analytically and compared to full 3D+1 pulse 
nonlinear propagation simulations. In the second section, 
experimental results are presented and compared to the- 
oretical results in several gases, leading to the determi- 
nation of the effective nonlinearity driving the ionization 
process in these gases. 



II. THEORETICAL APPROACH 

The method presented in this paper is based on the 
fact that any nonlinear process is enhanced when, at 
constant energy level, two coherent pulses spatially in- 
terfere as compared to a single pulse experiment. In the 
first subsection, the principle is demonstrated with in- 
finite plane waves. In the second subsection, the con- 
cept is extended analytically to gaussian beams propa- 
gating linearly. Finally, the result is compared to real- 
istic full 3D+1 nonlinear propagation simulations. The 
comparison between analytical calculations and numeri- 
cal simulations is of particular importance when dealing 
with bulk media since nonlinear effects usually lead to 
strong deviations between analytical results obtained in 
the linear regime and realistic propagation resulting in 
large uncertainty on the experimental determination of 
the nonlinearity. 



A. Plane waves case 

In this section, analytical gain on a nonlinear process 
is derived when two crossing plane waves with same car- 
rier frequencies interfere. Assuming infinitely extended 
plane waves, the electric field envelope e composed of 
two distinct crossing fields writes 
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where £i (resp. £2) is the electric field envelope of the 
first (resp. second) pulse and (j) is the relative phase be- 
tween the two pulses. If one defines e\ (resp. 82) such 
that l^ip = Ii (resp. \£i\^ = Ii) where Ii (resp. I2) is 
the intensity of the first (resp. second) beam, the total 
intensity then writes 



/ = /l + /2 + 2£i£2 COS (j) 

= It (^l -\- 2y/a{l — a) cos( 
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with It = /i +^2 the total intensity assuming no interfer- 
ence and a = Ii/It the relative intensity of the first pulse 
as compared to It- Assuming that the studied nonlinear 
process V is directly proportional to the N^^ power of 
intensity, the averaged signal over an interference fringe 
writes 
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with Ca = 2^ya{l — a) the fringe contrast. The signal 
Sl induced by a single beam with the same averaged 
intensity is given by 



The gain Ga then writes 
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where ^Fi is the hypergeometric function. Figure [TJa) 
shows the gain Ga as a function of a for several non- 
linearities. The gain on the nonlinear process can easily 
be controlled by adjusting the relative energy of the two 
pump beams. The maximal gain is obtained when the en- 
ergy is equally shared between the two beams {a = 0.5), 
i.e. when the fringe contrast is maximal (C = 1). In 
this case, the gain G1/2 compared to a single pulse of 
intensity It writes 
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where T : z t^~^ dt denotes the gamma func- 

tion. As shown in Fig. [Hb), the obtained gain is mono- 
tonic with the nonlinearity so that the determination 
of the gain leads in turn to the unambiguous determina- 
tion of the nonlinearity. Using this method, the deter- 
mination of the nonlinearity only requires two measure- 
ments {Sl and S1/2) unlike conventional methods which 
need to reconstruct the full dependence of the nonlinear 
process as a function of the intensity. 




Figure 1. (Color online) (a) Gain on the nonlinear process as 
a function of a for different nonlinearity orders N. Optimal 
gain on the nonlinear process (b) and energy gain (c) as a 
function of the nonlinearity. 



Note that an equivalent nonlinear signal would be gen- 
erated with a single pulse of intensity Ieq=G^lT where 
Ge is the energy gain. One can show that 
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As a consequence, the overall energy needed for the 
production of a nonlinear process in the two-beam geom- 
etry can be reduced up to a factor 2 as compared to the 
single beam geometry. 



B. Gaussian waves case 

In the previous section, the gain on the generation of 
any nonlinear process has been derived for infinite plane 
waves. Considering plane waves is a useful approxima- 
tion for understanding the proposed method and pro- 
vides very simple analytical formula. Nevertheless, most 
of nonlinear processes are produced with non collimated 
beams. In this section, the gain on a nonlinear process is 
therefore derived in the case of gaussian beams undergo- 
ing a linear propagation. The results are then compared 
to the gain numerically calculated with a 3D+1 nonlinear 
propagation equation. 



1. Analytical considerations 

In the case of two gaussian beams of wavenumber /cq 
crossing with an angle 2^, the intensity profile in the 
bisecting plane (x,z) as a function of the propagation 
distance z is given by 

I{X, Z) = /i(x, Z) + I2{X, Z) + /cross (^, z) COS (2/Co siufe), 

(8) 

with 
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A{z) = 
Rz = 
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where Aq is the wavelength, qz is the complex radius, ctq 
is the beam waist, /t is the total intensity, and Zr is the 
Rayleigh length. In the particular case a=0.5, it reduces 
to 



Figure 2. (Color online) (a) Gain on the nonlinear process as 
a function of a for several nonlinearity orders in the case of 
gaussian beams. Optimal gain on the nonlinear process (b) 
and energy gain (c) as a function of the nonlinearity. 



(about one order of magnitude for nonlinearity A/'=8-9) 
of the nonlinear process generation, in particular if one 
compares with standard control strategy based on optical 
pulse shaping Q. 
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with B = ^^2°/^? — ^^0 sin 0. 

The gain on the nonlinear process integrated over the 
propagation distance is given by 



fjl^{x,z)dxdz 
^11 {x^ z)dxdz 
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Figure [2fa) shows the gain as a function of a for 
several nonlinear it ies. The initial conditions were cho- 
sen to be Ao=800 nm, (Jq=22 /im, and 6>=3°. The gain 
integrated over the propagation length z is considerably 
reduced as compared with the plane wave analytical solu- 
tions. This reduction is diminished when the interaction 
region is limited by the length of the nonlinear medium 
itself rather than by the crossing effect. For low nonlin- 
earities, the generation of the nonlinear process is even 
reduced as compared with the single beam case. Figures 
[2jb,c) display G1/2 and the energy gain Ge as a function 
of the nonlinearity order calculated with optimal initial 
conditions, i.e. a = 0.5. As it is the case for plane 
waves, one can show that Ge 
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2. As a consequence, 

the same amount of nonlinear signal can be obtained in 
the two-pulse scheme with a total energy down to half 
the one required in the single pulse geometry. Moreover, 
the presented method leads to very efficient optimization 



2. Numerical considerations 



We chose to test the presented method with the pho- 
toionization process. In this case, the interferences be- 
tween the two pulses lead to the formation of a plasma 
grating. Such a grating has been the subject of rising 
interest. For instance, it has been shown that a mov- 
ing plasma grating formed by two intersecting lasers of 
slightly different central frequency leads to an energy ex- 
change between them [5| or that a stationary grating 
formed with pulses of same frequency results in the di- 
vergence or in the redirection of a third laser beam [6[. 
In order to assess the analytical results derived above 
and evaluate the influence of nonlinear propagation ef- 
fects such as intensity clamping in realistic experimental 
conditions, full 3D+1 numerical propagation calculations 
were performed. The propagation in both single (a = 1, 
£^ = 50 /iJ) and two-beam (a = 0.5, Ei = E2 = 25 /iJ) 
geometries were considered, allowing the comparison of 
the amount of electrons generated ah along the propaga- 
tion. The equation driving the propagation of the global 
electric field envelope e reads [8] 



dzS = i{kz- kiuj) £^ ^ (^Pnl 



2eoc' 



-J 



-^losses 5 

(13) 



4 



with 



k{co) = n{co)oj/c, = y/fc2H-(fc2+fc2), (14) 
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where u is the angular frequency, n is the refractive index, 
c is the hght velocity, rUe and e are the electron mass and 
charge respectively, n2 is the nonlinear refractive index, 
Pat is the atoms density, N is the ionization nonlinearity, 
cr AT is the ionization cross-section, and is the effective 
collision frequency. We chose N=7.5 since it is the widely 
admitted effective nonlinearity driving the ionization rate 
in argon. Finally, we considered the same gas pressure 
(0.4 bar) than the one used during our experiments. 

In the two-beam case, two gaussian pulses propagating 
slightly off- axis {0 = ±3°) with respect to the z axis 
and crossing along the x dimension have been considered. 
The initial electric field then writes as 



e{x,y,t,z = 0) = F{t)G{y){e+H+{x)+e-H-{x)) (15) 
with 
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and (70=22 pm, 6^=3°, Ao=800 nm, ^o=-4 mm, and (7^=85 
fs, and E±=25 jaJ. 



all along the propagation as shown in Figure [3l^b). The 
gain of free electrons averaged over the whole propagation 
obtained analytically (Eq. [T2]) and numerically only dif- 
fers by 0.2 % (5.45 and 5.44, respectively). It ensures that 
nonlinear propagation effects do not significantly alter 
the linear propagation, at least in the present initial en- 
ergy and pressure conditions, and consequently validates 
the analytical results. The proposed method is there- 
fore well suited for fast and accurate determinations of 
unknown nonlinearity for experimental conditions where 
nonlinear propagation effects remains negligible. 
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Figure 3. (Color online) (a) Amount of free electrons per 
unit length as a function of the propagation distance calcu- 
lated according to the full 3D propagation equation in the one 
beam (dashed blue) and two-beam configurations (solid red), 
(b) Comparison between analytical and numerical gain as a 
function of the propagation distance. 



III. EXPERIMENTAL RESULTS 



S. Discussions 

Figure [3l^a) displays the amount of free electrons per 
length unit generated in the single pulse case (blue 
dashed line) and in the two-pulse configuration (solid red 
line) calculated with the help of the numerical simula- 
tions. As expected, the production of free electrons is 
enhanced as soon as the pulses overlap. The maximal 
gain, which is obtained when the two pulses perfectly 
overlap (^=0), is about 37 in perfect agreement with Eq. 
[6l which indicates that the plane wave approximation is 
valid in this particular situation. A remarkable agree- 
ment between analytical and numerical results is found 
when comparing the evolution of the free electrons gains 



A. Experimental setup 

The setup is depicted in Fig. HI The optical source is 
a 1 kHz amplified femtosecond laser delivering vertically 
polarized 600 /iJ, 100 fs pulses at 796.3 nm. The beam is 
split into two parts, the relative delay r between the two 
paths (Pi and P2) being controlled by means of a mo- 
torized delay line. The energy in the two optical paths 
is adjusted by means of an half- wave plate and a polar- 
izer. In the path P2, a Mach-Zehnder interferometer is 
inserted to produce two co-propagating pulses. The first 
one is used as a pump during the two-pulse experiments 
and is blocked for single pump experiments while the sec- 
ond one, delayed by around ri=500 fs, is used to probe 
the plasma density generated by the single or the two- 
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pump beam, depending on the experiment. The energy 
of the probe is about 4 % of the total energy in P2. More- 
over, its polarization is rotated by 90° with respect to the 
polarization of the pumps. The three pulses are then fo- 
cused in a static cell. Pi and P2 are crossed with an 
angle 26 ~5.2°. The waist of the pulses, aexp = 22 /im, 
were measured at weak intensity with the help of a cam- 
era. After the static cell, the probe is selected with the 
help of a polarizer. The plasma density is measured by 
the cross-defocusing technique as described in [3, [9| . As 
shown in 0, Q in the case of parabolic pulses, the de- 
focusing signal is proportional to An^, i.e. the peak 
to valley change of refractive index experienced by the 
probe beam. This is confirmed in the gaussian pulse case 
by our full 3D+1 numerical calculations simulating the 
experimental defocusing setup. In this experiment, the 
cross-defocusing signal is then proportional to the square 
of refractive index change resulting from the ionization 
mechanism and accordingly proportional to the square of 
the amount of free electrons generated during the exper- 
iment. It then allows a direct experimental measurement 
of the latter. Since the defocusing signal induced in the 
two-pulse scheme is about 30 times higher than the one 
induced during the single pulse experiment, a calibrated 
optical density is added before the photomultiplier dur- 
ing the former in order to keep the detection dynamics 
constant and to avoid a saturation of the photomultiplier. 



line) and two-pulse (25 fiJ-\-25 /iJ, resp. 40 /iJ-h40 /iJ, 
red line) cases as a function the pump-probe delay. Note 
that it was checked with the help of 3D+1 numerical 
propagation simulations that intensity clamping does not 
occur at this pressure up to 90 TW/cm^. As far as the 
single pulse case is concerned, the defocusing signal in- 
creases as soon as the pump and probe pulses temporally 
overlap and then remains roughly constant for the rest 
of the scan, as expected from plasma-induced defocusing 
(both electrons recombination and diffusion are negligible 
over the considered temporal window). In the two-pulse 
experiment, the defocusing signal remains hardly distin- 
guishable on the scale of the figure until the two pulses 
temporally overlap, resulting in a sharp increase of the 
signal. The gain on the production of free electrons is 
retrieved from the experimental curves by comparing the 
post-pulse defocusing signal induced by a single pump to 
the one induced when the two pumps perfectly tempo- 
rally overlap. The result obtained at £^tot = 50/iJ as a 
function of a is depicted in Fig. [5fc). A nonlinearity 
A^eff=7.4± 0.05 is found when fitting the experimental 
data with the analytical model (Eq. [T2]) . in good agree- 
ment with previous works performed in the same inten- 
sity regime. 




Figure 4. (Color online) Top view of the experimental setup. 
BS: beam splitter. GP: Glan cube polarizer. OD: optical 
density Pi (P2): path 1 (2). S: probe. 

Note that, in the two-beam case, the defocusing sig- 
nal in the probe propagation direction is sensitive to the 
plasma density averaged over an interference fringe (as 
calculated in the analytical section) while the modula- 
tion due to interferences is responsible for the redirection 
of a small part of the probe 



B. Results and discussion 

Figure [5ja), resp. (b), shows the defocusing signal in 
argon (0.4 bar) in both one (50 /iJ, resp. 80 /iJ, black 
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Figure 5. (Color online) Defocusing signal in argon at 50 
/iJ (a) and at 80 /x J (b) as a function of the delay between 
pump and probe beams for the single and two-beam pump 
experiments, (c) Free electrons gain as a function of a in 
argon at 50 /iJ. Experiments 1, 2, and 3 correspond to three 
independent measurements sets. Measured nonlinearity as a 
function of laser peak intensity for (d) argon, (e) ethylene, 
and (f) krypton. 

The above effective nonlinearity is intensity dependent. 
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When an electron oscillates in the field, it acquires a pon- 
deromotive energy that prevents the electron to escape 
from the atomic potential, making the ionization harder. 
Furthermore, as the electric field increases tunnel ion- 
ization takes over resulting in a deviation from the mul- 
tiphoton ionization regime. Both effects are expected 
to reduce the effective nonlinearity as the intensity in- 
creases. In order to study such a dependence, the global 
energy was changed from about 20 /iJ up to 80 fiJ. Fig- 
ures[5Kd-f) show the nonlinearity of the ionization process 
as a function of the global energy in argon, ethylene, and 
krypton. These gases have been chosen because their 
ionization potential /p differs by several photon energy 
units (see Table IH) (Eph ^ 1.56 eV). For ah gases, the 
effective nonlinearity decreases by more than one unit 
as the intensity increases from 20 to 80 TW/cm^ and 
can be fitted by a linear function A^exp(^) = + 
with respect to intensity /. The coefficients of the fit 
are summarized in Table HI Finally, one has to empha- 
size that the nonlinearity measurements are performed 
by two measurements only (one and two-pulse configura- 
tion) unlike conventional technics. This technique then 
provides a very convenient method to study highly non- 
linear systems. 





C2H4 


Kr 


Ar 


Ip (eV) 


10.51 


14.00 


15.76 


Nth@796.3 nm 


7 


9/10^*^ 


11 


Nexp @ 50 TW/cm^ 


5.7 


7.5 


7.4 


lo (TW/cm^) 


27.3 


26.9 


37.3 




7.6 


9.3 


8.7 



Table I. Comparison between the theoretical nonlinearity in 
the pure multiphoton regime and the nonlinearity experimen- 
tally measured. *^*'* In the case of Kr, the theoretical nonlin- 
earity at 796.3 nm is 9 while it is 10 at 797 nm, both wave- 
lengths being covered by the laser bandwidth. 

These results could be used to improve the numerical 



models of ionization used, for instance, in filamentation 
studies. Nevertheless, since the presented method only 
gives the nonlinearity of the process but is unable to de- 
termine the absolute ionization probability, it should be 
coupled to a self-referenced method. For instance, by 
comparing the ionization induced refractive index change 
with the one induced by molecular alignment [4] , it could 
be possible to determine the absolute cross-section of the 
ionization process. 

IV. CONCLUSION 

In conclusion, a simple and fast all-optical method for 
measuring nonlinearities and optimizing their concomi- 
tant process is presented. The principle of the present 
method is validated both analytically and numerically. 
It is applied to the nonlinear photoionization of several 
gases. Moreover, it is shown that interferences between 
two crossing pulses lead to a net gain that can reach 
several order of magnitude on the nonlinearity genera- 
tion efficiency at constant energy level. Finally, one has 
to emphasize that larger enhancement can be achieved 
by using more that two pulses. For instance, the maxi- 
mal gain expected in the three plane waves configuration 
is ^3 beams = 3"^ (1/2 - AT, -AT, 1, 4), the energy 
needed to produce a nonlinear effect being reduced up to 
a factor three. This improvement would require a stabi- 
lized interferometer for maximizing multi- waves interfer- 
ences. 
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